Scale Invariance and Lack of Self- Averaging in Fragmentation 
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We derive exact statistical properties of a class of recursive fragmentation processes. We show 
that introducing a fragmentation probability < p < 1 leads to a purely algebraic size distribution 
in one dimension, P{x) oc x~^^. In d dimensions, the volume distribution diverges algebraically 
in the small fragment limit, P{V) ~ V ' with 7 — 2p^''*. Hence, the entire range of exponents 
allowed by mass conservation is realized. We demonstrate that this fragmentation process is non- 
self- averaging. Specifically, the moments Y^ — ^ - a;f exhibit significant fluctuations even in the 
thermodynamic limit. 
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Numerous physical phenomena are characterized by a 
set of variables, say {xj}, which evolve according to a 
random process, and are subject to the conservation law 
^ ■ Xj — const. An important example is fragmentation, 
with applications ranging from geology |l|] and fracture 
H to the breakup of liquid droplets Q and atomic nu- 
clei [IIJ^]. Other examples include spin glasses [g|, where 
Xj represents the equilibrium probability of finding the 
system in the j*** valley, genetic populations, where Xj is 
the frequency of the j*'' allele HM] , and random Boolean 
networks [|jlO|. 

In most cases, stochasticity governs both the way in 
which fragments are produced, and the number of frag- 
mentation events they experience. For example, in frag- 
mentation of solid objects due to impact with a hard 
surface fragments may bounce several times before com- 
ing to a rest pd| . The typical number of fragmenta- 
tion events may vary greatly as it depends on the ini- 
tial kinetic energy. Another seemingly unrelated exam- 
ple is DNA segmentation algorithms where homogeneous 
subsequences are produced recursively from a inhomoge- 
neous sequence until a predefined homogeneity level is 
reached [p2| . Here, the number of segmentation events is 
determined by the degree of homogeneity of the original 
sequence. 

In this study, we examine fragmentation models with 
two types of objects - stable towards fragmentation and 
unstable. We show that the size distribution is algebraic, 
and that the entire range of power-laws allowed by the 
underlying conservation laws can be realized by tuning 
the fragmentation probability. Additionally, such pro- 
cesses are characterized by large sample to sample fluc- 
tuations, as seen from analysis of the moments of the 
fragment size distribution. 

Specifically, we consider the following recursive frag- 
mentation process. We start with the unit interval and 
choose a break point Hn [0, 1] with a uniform probability 
density. Then, with probability p, the interval is divided 
into two fragments of lengths I and l~l, while with prob- 
ability q = 1 —p, the interval becomes "frozen" and never 
fragmented again. If the interval is fragmented, we recur- 



sively apply the above fragmentation procedure to both 
of the resulting fragments. 

First, let us examine the average total number of frag- 
ments, N. With probability q a single fragment is pro- 
duced, and with probability p the process is repeated 
with two fragments. Hence N — q + 2pN, yielding 



(q/il-2p), ifp<l/2; 
\oo, ifp>l/2. 



(1) 



The average total number of fragments becomes infinite 
at the critical point pc = 1/2, refiecting the critical na- 
ture of the underlying branching process ||l^ . 

Next, we study P{x), the density of fragments of length 
X. The recursive nature of the process can be used to ob- 
tain the fragment length density 



P{x) = q6{x -l) + 2p f ^p(- 

Jx y \y 



(2) 



The second term indicates that a fragment can be created 
only from a larger fragment, and the y^^ kernel reflects 
the uniform fragmentation density. Eq. (0) can be solved 
by introducing the Mellin transform 



M{s) = / dxx'~^P{x). 



(3) 



Eqs. (|) and (|) yield M{s) =q + 2ps-^M{s) and as a 
result 



M{s) = q 



2p 



s~2p 



(4) 



The average total number M(l) = N is consistent with 
Eq. (|l|), and the total fragment length M(2) = 1 is con- 
served in accord with 1 = J dx xP{x). (Here and in the 
following the integration is carried over the unit interval, 
i.e., < a; < 1.) The inverse Mellin transform of Eq. (0) 
gives 



P(x) ^q[S{x- l) + 2p: 



-2pl 



(5) 



Apart from the obvious (5-function, the length density is 
a purely algebraic function. In particular, the fragment 
distribution diverges algebraically in the limit of small 
fragments. Given such an algebraic divergence near the 
origin P{x) ^ x^'^ , length conservation restricts the ex- 
ponent range to 7 < 2. In our case 7 = 2p, and since 
< p < 1, the entire range of acceptable exponents 
emerges by tuning the only control parameter p. 

Interestingly, at the critical point Pc = \^ the frag- 
ment length distribution becomes independent of the ini- 
tial interval length. Starting from an interval of length 
L, Eq. (ph can be generalized to yield 



P{x) = q5{x -L) + 2pqL^-^Px-'^P. 



(6) 



Thus, the critical point may be detected by observing 
the point at which the segment distribution becomes in- 
dependent of the original interval length. 

The recursive fragmentation process can be general- 
ized to d dimensions. For instance, in two dimensions we 
start with the unit square, choose a point {xi, X2) with a 
uniform probability density, and divide, with probability 
p, the original square into four rectangles of sizes xi XX2, 
xi X (1 -X2), (l-xi) X X2, and (l-xi) x (l-a;2). With 
probability q, the square becomes frozen and we never 
again attempt to fragment it. The process is repeated 
recursively whenever a new fragment is produced. 

Let ^(x), X = (xi, . . . ,Xd), be the probability density 
of fragments of size xi x ■ ■ ■ x Xd- This quantity satisfies 



P(x) = q6{x - 1) 



2'^p 



yi---yd V 2/1 ' " ' ' yd 



(7) 



with J dy — J dyi ■ ■ ■ J dya- Following the steps leading 
to Eq. (0), we find that the d-dimensional Mellin trans- 
form, defined by Af (s) = / rfx^i'"^ • • • x^''"^P(x) with 
the shorthand notation s = (si, . . . , Sd) obeys 



M(s) 



si ■ 



■ Srf - a" 



with a = 2p^''^. (8) 



Eq. (g) gives the total average number of fragments, 
N = Af (1) = q/{l ~ 2'^p) a p < 2-"^ and A^ = 00 if 
p > 2^'^. One can also verify that the total volume 
A/(2) = 1 is conserved. Interestingly, there is an ad- 
ditional infinite set of conserved quantities: all moments 
whose indices belong to the hyper-surface s| • ■ • sjj = 2'' 
satisfy A/(s*) = 1. In a continuous time formulation of 
this process the same moments were found to be integrals 
of motion pJ-M . The existence of an infinite number of 
conservation laws is surprising, because only the volume 
conservation has a clear physical justification. 

Next, we study the volume density P{V)^ defined by 



P{V)= fd^P{^)5{V-xi---Xd) 



(9) 



The Mellin transform M{s) = J dVV''-'^P{V) can be ob- 
tained from Eq. (Sj by setting Si — s, 



M{s) = q 



(10) 



Using the d'^ root of unity, C = e^'^'/^, and the iden- 



tity 



I Y^d-l 

d Z^fc=o 



1 _c_ 



M(s) can be expressed as a 



s<'-l d Z.^fc=o s-C'' ' 

sum over simple poles at aC,^ . Consequently, the inverse 
Mellin transform is given by a linear combination of d 
power laws 



P{V) = q 



d-l 



s{v-i) + ^Y.("y''''^ 



fc=0 



(11) 



One can verify that this expression equals its complex 
conjugate and hence, it is real. Additionally, the one- 
dimensional case (||) is recovered by setting d — 1. 

The small-volume tail of the distribution can be ob- 
tained by noting that the sum in Eq. (O) is dominated 
by the first term in the series, which leads to 



P{V) ~ AV-'' 



V^O, 



(12) 



with ^ = a — 2p^l'^ and A ~ aq/d. Although the value 
of the exponent changes, the possible range of exponents 
for this process remains the same since < 2p^/''- < 2 
when < p < 1. In the infinite dimension limit, P(V) 
becomes universal: P{V) ~ V~'^. 

The leading behavior of P{V) in the large size limit 
can be derived by using the Taylor expansion and the 

identity J2k=o C^" ~ '^".0 for n = 0, . . . , d — 1. One finds 
that in higher dimensions the volume distribution van- 
ishes algebraically near its maximum value, 



p{v)^Bd{i-vy 



V 



(13) 



withBd = a'^/(d-l)!. 

In fact, the entire multivariate fragment length den- 
sity can be obtained explicitly. This can be achieved by 
expanding the geometric series 



d / \ n-l-1 



si---Sd 



n>Oi=l 



and performing the inverse Mellin transform for 
each variable separately. Using the transform 

J dxx^^^ [l^i] — nls~"~^ gives 

P(x)=q[J(x-l) + a'^Fd(z)], (14) 

with the shorthand notations 

l/d 



^*>-i:^ 



n=0 



ni4 



(15) 



In two dimensions, ^2(2) = Iq{2z) where /q is the modi- 
fied Bessel function. 



The small size behavior of -P(x) can be obtained by 
using the steepest decent method. The leading tail 
behavior, Fd{z) ~ (27rz)T~e^'' for z ^ 1, corresponds 
to the case when at least one of the lengths is small, 
i. e. Xi <C 1. Returning to the original variables we see 
that the fragment distribution exhibits an unusual "log- 
stretched-exponential" behavior 



P(x) 



d 

ni4 



exp 




(16) 



The fragment distribution represents an average over 
infinitely many realizations of the fragmentation process; 
hence, it does not capture sample to sample fluctuations. 
These fluctuations are important in non-self-averaging 
systems, where they do not vanish in the thermodynamic 
limit. Useful quantities for characterizing such fluctua- 
tions are the moments Yq ]l7|,P^ 



Yr, 



T.- 



(17) 



when a, f3,a + (3 > 2p — 1, and (YaYp) — oo otherwise. 

Eq. (|l|) shows that (YaYp) 4- {Ya){Yfi), and in par- 
ticular, (Y^) ^ (Ya)"^ . Therefore, fluctuations in Yq do 
not vanish in the thermodynamic limit, and the recursive 
fragmentation process is non-self-averaging. While for 
p < 1/2 non-self-averaging behavior is expected because 
the average number of fragments is finite, the emergence 
of non-self-averaging for p > 1/2 is surprising. Hence, 
statistical properties obtained by averaging over all real- 
izations are insufficient to probe sample to sample fluc- 
tuations. 

In principle, higher order averages such as (Y^) can be 
calculated recursively by the procedure outlined above. 
The resulting expressions are cumbersome and not terri- 
bly illuminating. Instead, one may study the distribution 
Qa{Y) which obeys 



Qc.{Y) ^ qS{Y - 1) 



(22) 



Y -Z 



(1 -/)""" V(i-0" 



:Q 



where the sum runs over all fragments. 

We are interested in the average values (Ya) and 
{YoYp). For integer a, (Yq) is the probability that a 
points randomly chosen in the unit interval belong to the 
same fragment. The expected value of Yq satisfies 



{Y^)=q+p{Y^) / dy [y" + (1 - y)"] 



(18) 



The first term corresponds to the case where the unit in- 
terval is not fragmented, and the second term describes 
the situation when at least one fragmentation event has 
occurred. Eq. (nsh gives 



{Yc) = q 



1 



2p 



a+l~2p 



(19) 



> 'Zp — i, and (Ya) = oo if a < 2p — 1. As expected, 
Eq. (^^ agrees with the moments of P{x) obtained by 
integrating Eq. (||), (Ya) = Jdxx°'P{x). 

However, higher order averages such as (YaY/s) do not 
follow directly from the fragment density. For integer a 
and /3, (YaYp) is the probability that, ii a + (3 points 
are chosen at random, the first a points all lie on the 
same fragment, and the last /3 points all lie on another 
(possibly the same) fragment. This quantity satisfies 



{YaYp} = q + p{YaYp) / dy [2/"+^ + {l-y) 



\a+fJ] 



(20) 



+ p{Ya){Yp) I dy [y"(l - yf + (1 - y)"/] , 
yielding 

(YaYp) = q + 

+ 2p 



2pq 

a + f3+l-2p 

r(a + i)r(/3 + i) (r„)(Y» 



(21) 



r(a + /3 + l) a + (3+l-2p 



In addition to the recursive nature of the process, we 
have employed extensivity, i.e, (Ya) oc L" in an interval 
of length L. 

Clearly, Yq = -^ and Yi = 1, and therefore Qi{Y) = 
6{Y — 1) and Qo{N) can be determined analytically as 
well. Generally, different behaviors emerge for a > 1 and 
a < 1. We concentrate on the former case where the sup- 
port of the distribution Qa{Y) is the interval [0,1]. The 
Laplace transform, RaW = Jq dY e^^^Qa{Y), obeys 

Ra {X)=qe-^+p f dl Ra [A/"] Ra [A(l -/)"]. (23) 
Jo 

The behavior of Qa{Y) in the limit F ^ is reflected 
by the asymptotics of i?a(A) as A ^ 00. Substituting 
_Rq,(A) ~ exp{—AX^) into both sides of Eq. (p3|), evalu- 
ating the integral using steepest decent and equating the 
left and right hand sides gives (3 = \/a. Consequently, 
we flnd that the distribution has an essential singularity 
near the origin 



}a{Ya) ~ exp 



-BYa 



Ya 



0. 



(24) 



Extremal properties can be viewed as an additional 
probe of sample to sample fluctuations. We thus con- 
sider C{x), the length density of the largest fragment. 
For a self-averaging process with an inflnite number of 
fragments, one expects C{x) — > 5{x) in the thermody- 
namic limit. To see that C{x) is non-trivial for any p, let 
us first determine C{x) for x > 1/2. In this region. 



C{x) = qS(x — 1) + p 



'dy 



(25) 



If the original unit interval has not been fragmented, the 
largest fragment is obviously the unit interval. If the first 
fragmentation is performed, only one of the two resulting 
fragments can be larger than x > 1/2. Therefore, only 
subsequent breaking of this fragment (of length y > x) 
can contribute to C{x), which explains Eq. (|25|). 

Eq. ( p5| ) is similar to Eq. (||) , and can be solved by the 
same technique to give 

C{x) =q6{x-l)+pqx-P for x > 1/2. (26) 

In the complementary case of a; < 1/2, C{x) satisfies 



.yj J y \y 

l-x 1/2 



'^<^'-/7<;)-/ 



d-y ^f ^ \ £ f^ 
y \^-yJ ^ \y 

1/2 



(27) 



The first term on the right-hand side of this equation is 
constructed as in Eq. (^5|) : if we first break the unit in- 
terval into two fragments of lengths y > 1/2 and 1 — y, 
then for 1 — y < a; the longest fragment is produced by 
breaking the fragment of length y. The next two terms 
describe the situation when I ~ y > x, so the longest 
fragment can arise out of breaking any of the two frag- 
ments. The factors C-{u) = J^ dv C{v) guarantee that 
the longest fragment of length x comes from the fragment 
of length V in the first generation. Since C{x) obeys dif- 
ferent equations in different regions, it looses analyticity 
on the boundaries. Namely, C{x) possesses an infinite set 
of singularities at x = 1/k which become weaker as k in- 
creases. Similar singularities undcrly extremal properties 
of a number of random processes including random walks, 
spin glasses, random maps, and random trees ]7|,p|,p7|-p^ . 

In summary, we have found that recursive fragmenta- 
tion is scale free, i.e., the fragment length distribution is 
purely algebraic. In higher dimensions, the volume distri- 
bution is a linear combination of d power laws, and conse- 
quently, an algebraic divergence characterizes the small- 
fragment tail of the distribution. A number of recent im- 
pact fragmentation experiments reported algebraic mass 
distributions with the corresponding exponents ranging 
from 1 to 2 [Q . It will be interesting to further examine 
whether our simplified model is appropriate for describ- 
ing fragmentation of solid objects. 

We have also found that the recursive fragmenta- 
tion process exhibits a number of features that arise in 
other complex and disordered systems, such as non-self- 
averaging behavior and the existence of an infinite num- 
ber of singularities in the distribution of the largest frag- 
ment. These features indicate that even in the thermo- 
dynamic limit sample to sample fluctuations remain, and 
that knowledge of first order averages may not be suffi- 
cient for characterizing the system. Our ID model is 



equivalent to applying the aforementioned DNA segmen- 
tation algorithm to a random sequence. It will be inter- 
esting to study self-averaging and extremal properties of 
DNA sequences, which are known to have commonalities 
with disordered systems. Indeed, if these subtle features 
are found for genetic sequences as well, this would sug- 
gest that much caution should be exercised in statistical 
analysis of DNA. 
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